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Abstract 

In this paper we show how non-radial motions, originating from 
the tidal interaction of the irregular mass distribution within and 
around protoclusters, can solve some of the problems of the CDM 
model. Firstly the discrepancy between the CDM predicted two- 
points correlation function of clusters and the observed one. We 
compare the two-points correlation function, that we obtain taking 
account of non-radial motions, with that obtained by Sutherland 
& Efstathiou (1991) from the analysis of Geller & Hucra's (1988) 
deep redshift survey and with the data points for the APM clusters 
obtained by Efstathiou et al. (1992). Secondly the problem of the 
X-ray clusters abundance over-production predicted by the CDM 
model. In this case we compare the X-ray temperature distribu- 
tion function, calculated using Press-Schechter theory and Evrard's 
(1990) prescriptions for the mass-temperature relation, taking also 
account of the non-radial motions, with Henry & Arnaud (1991) and 
Edge et al. (1990) X-ray temperature distributions for local clus- 
ters. We find that in both cases the model is in good agreement with 
experimental data. Finally we calculate the bias coefficient using a 
selection function that takes into account the effects of non-radial 
motions, and we show that the bias so obtained can account for a 
substantial part of the total bias required by observations on cluster 
scales. 

keywords: cosmology: theory-large scale structure of Universe - 
galaxies: formation 
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1 Introduction 



Although at his appearence the standard form of CDM was very success- 
ful in describing the observed structures in the universe (galaxy clustering 
statistics, structure formation epochs, peculiar velocity flows) (Peebles, 
1982; Blumenthal et al. 1984; Bardeen et al. 1986; White et al. 1987; 
Frenk et al. 1988; Efstathiou 1990) recent measurements have shown sev- 
eral deficiencies of the model, at least if any bias of the distribution of 
galaxies relative to the mass is constant with scale. Some of the most diffi- 
cult problems to reconcile with the theory are the strong clustering of rich 

clusters of galaxies, £ C c(?") — ( 25fe _i Mpe ) , far in excess of CDM predictions 
(Bahcall & Soneira 1983), the X-ray temperature distribution function of 
clusters, over-producing the observed cluster abundances (Bartlett & Silk 
1993), the conflict between the normalisation of the spectrum of the per- 
turbation which is required by different types of observations. 
Alternative models with more large-scale power than CDM have been in- 
troduced in order to solve the latter problem. Several authors (Peebles 
1984; Efstathiou, Sutherland & Maddox 1990, Turner 1991, White, Efs- 
tathiou & Frenk 1993) have lowered the matter density under the critical 
value (Q m < 1) and they have added a cosmological constant in order to 
retain a flat universe (Q m + Q\ = 1) . The spectrum of the matter density 
is specified by the transfer function, but its shape is affected because of the 
fact that the epoch of matter-radiation equality (characterized by a red- 
shift z eq ) is earlier, 1 + z eq being increased by a factor l/f2 m . Around the 
epoch z\, where z\ = (tt m /n^y/ 3 — 1, the growth of the density contrast 
slows down and ceases after z\. As a consequence the normalisation of the 
transfer function begins to fall, even if its shape is retained. Mixed dark 
matter models (MDM) (Bond et al. 1980; Shafi & Stecker 1984; Valdarnini 
& Bonometto 1985; Schaefer et al. 1989; Holtzman 1989; Schaefer 1991; 
Shaefer and Shafi 1993; Holtzman & Primack 1993) increase the large-scale 
power because neutrinos free-streaming damps the power on small scales. 
Alternatively changing the primeval spectrum several problems of CDM 
are solved (Cen et al. 1992). Finally it is possible to assume that the 
threshold for galaxy formation is not spatially invariant but weakly modu- 
lated (2% — 3% on scales r > lOh^Mpc) by large scale density fluctuations, 
with the result that the clustering on large-scale is significantly increased 
(Bower et al. 1993). 

Here we propose a different solution to several of the CDM model problems 
connected to the non-radial motions developing during the protocluster 
evolution. 

It has long been speculated that angular momentum could have a fun- 
damental role in determining the fate of collapsing proto-structures and 
several models have been proposed in which the galaxy type can be cor- 
related with the angular momentum per unit mass of the structure itself 
(Faber 1982; Kashlinsky 1982; Fall 1983). Some authors (see Barrow & 
Silk 1981, Szalay & Silk 1983 and Peebles 1990) have proposed that non- 
radial motions would be expected within a developing proto-cluster due to 
the tidal interaction of the irregular mass distribution around them, typi- 
cal of hierarchical clustering models, with the neighboring proto-clusters. 
The kinetic energy of this non-radial motions opposes the collapse of the 
proto-cluster, enabling the same to reach statistical equilibrium before the 



final collapse (the so called previrialization conjecture by Davis & Peebles 
1977, Peebles 1990). This effect may prevent the increase of the slope of the 
mass autocorrelation function at separations given by £(r, t) ~ 1, expected 
in the scaling solution for the growth of £(r, t) but not observed in the 
galaxy two-point correlation function. The role of non-radial motions has 
been pointed by several authors (see Davis & Peebles 1983, Gorski 1988, 
Groth et al. 1989, Mo et al. 1993, Weygaert & Babul 1994, Marzke et al. 
1995 and Antonuccio & Colafrancesco 1997). Antonuccio & Colafrancesco 
derived the conditional probability distribution f p k(y\v) of the peculiar ve- 
locity around a peak of a Gaussian density field and used the moments of 
the velocity distribution to study the velocity dispersion around the peak. 
They showed that regions of the proto-clusters at radii, r, greater than the 
filtering length, Rf, contain predominantly non-radial motions. 
Non-radial motions change the energetics of the collapse model by intro- 
ducing another potential energy term. In other words one expects that 
non-radial motions change the characteristics of the collapse and in partic- 
ular the turn around epoch, t m , and consequently the critical threshold, S c , 
for collapse. One expects that non-radial motions produce firstly a change 
in the turn around epoch, secondly a new functional form for S c , thirdly 
a change of the mass function calculable with the Press- Schechter (1974) 
formula and consequently of the predicted X-ray temperature distribution 
function of clusters and finally a modification of the two-point correlation 
function. Moreover this study of the role of non-radial motions in the col- 
lapse of density perturbations can help us to give a deeper insight on the 
so called problem of biasing. As pointed out by Davis et al. (1985), unbi- 
ased CDM suffers of several problems: pairwise velocity dispersion larger 
than the observed one, galaxy correlation function steeper than observed 
(see Liddle & Lyth 1993 and Strauss & Willick 1995). The remedy to these 
problems is the concept of biasing (Kaiser 1984), i.e. that galaxies are more 
strongly clustered than the mass distribution from which they originated. 
The physical origin of such biasing is not yet clear even if several mecha- 
nisms have been proposed (Rees 1985; Dekel & Rees 1987; Dekel & Silk 
1986; Carlberg 1991; Cen & Ostriker 1992; Bower et al. 1993; Silk & Wyse 
1993). Recently Colafrancesco, Antonuccio & Del Popolo (1995, hereafter 
CAD) have shown that dynamical friction delays the collapse of \ow-v 
peaks inducing a bias of dynamical nature. Because of dynamical friction 
under-dense regions in clusters (the clusters outskirts) accrete less mass 
with respect to that accreted in absence of this dissipative effect and as a 
consequence over-dense regions are biased toward higher mass (Antonuccio 
& Colafrancesco 1995 and Del Popolo & Gambera, 1996). Non-radial mo- 
tions acts in a similar fashion to dynamical friction: they delay the shell 
collapse consequently inducing a dynamical bias similar to that produced 
by dynamical friction. This dynamical bias can be evaluated defining a 
selection function similar to that given in CAD and using Bardeen, Bond, 
Szalay and Kaiser (1986, hereafter BBKS) prescriptions. The plan of the 
paper is the following: in §2 we obtain the total specific angular momen- 
tum acquired during expansion by a proto-cluster. In §3 we find the effect 
of non-radial motion on the critical density threshold, S c . In §4 we derive a 
selection function for the peaks giving rise to proto-structures while in §5 
we calculate some values for the bias parameter, using the selection func- 
tion derived, on three relevant filtering scales. In §6 we find the effects 



of non-radial motions on the X-ray temperature distribution function and 
then we compare this prevision to the X-ray observed data. In § 7 we 
study how non-radial motions influence the cluters two-points correlation 
function. §8 is devoted to conclusions and discussions. 



2 Tidal torques in clusters evolution. 

The explanation of galaxies spins gain through tidal torques was pioneered 
by Hoyle (1949) in the context of a collapsing protogalaxy. Peebles (1969) 
considered the process in the context of an expanding world model show- 
ing that the angular momentum gained by the matter in a random co- 
moving Eulerian sphere grows at second order in proportion to t 5 ^ 3 (in a 
Einstein-de Sitter universe), since when the proto-galaxy is still a small 
perturbation, while in the nonlinear stage the growth rate of an oblate 
homogeneous spheroid decreases with time as t' 1 . More recent analytic 
computations (White 1984, Hoffman 1986, Ryden 1988) and numerical sim- 
ulations (Barnes & Efstathiou 1987) have re-investigated the role of tidal 
torques in originating galaxies angular momentum. In particular White 
(1984) expanded an analysis by Doroshkevich (1970) showing that the an- 
gular momentum of galaxies grows to first order in proportion to t and that 
the result of Peebles is a consequence of the spherical symmetry imposed 
in the model. White showed that the angular momentum of a Lagrangian 
sphere does not grow either in first or in second order while the angular 
momentum of a non-spherical volume grows to first order in agreement to 
Doroshkevich's result. Hoffman (1986) has been much more involved in 
the analysis of the correlation of the growth of angular momentum with 
the density perturbation S(r). He found an angular momentum- density 
anticorrelation: high density peaks acquire less angular momentum than 
low density peaks. One way to study the variation of angular momen- 
tum with radius in a galaxy is that followed by Ryden (1988). In this 
approach the protogalaxy is divided into a series of mass shells and the 
torque on each mass shell is computed separately. The density profile of 
each proto-structure is approximated by the superposition of a spherical 
profile, S(r), and a random CDM distribution, e(r), which provides the 
quadrupole moment of the protogalaxy. To first order, the initial density 
can be represented by: 

p(r)=p b [l + 6(r)][l + e(r)] (1) 
where pb is the background density and e(r) is given by: 

<M 2 > = P(k) (2) 

being P(k) the power spectrum, while the density profile is (Ryden & Gunn 
1987): 
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where v is the height of a density peak, £(r) is the two points correlation 
function, 7 and are two spectral parameters (BBKS, Eq. 4.6a, 4.6d) 



while $(7^, 7) is a function given in BBKS (Eq. 6.14). As shown by Ryden 
(1988) the net rms torque on a mass shell centered on the origin of internal 
radius r and thickness 5r is given by: 

<M 2 > 1/2 = VSO (yG) [(a 2m (r) 2 )(g 2m (r) 2 > - (a 2m (r)g 2 * m (r)) 2 ] 1/2 (4) 

where g; m , the multipole moments of the shell and ai m , the tidal moments, 
are given by: 
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(a 2m (r) 2 ) = /" dkP (A;) ^ (fc r ) 2 ( 6 ) 



r 



(« 2m (r)g 2 * m (r)) = ^jP 6 M sft j kdkP (k) 3l (kr) j 2 (kr) (7) 

where M s /j is the mass of the shell, j%(r) and j 2 (r) are the spherical Bessel 
function of first and second order while the power spectrum P{k) is given 
by: 

P(k) = Ak' 1 [ln(l + 4.164A;)] 2 
(l92.9 + 1340A; + 1.599 x 10 5 k 2 + 1.78 x 10 5 k 3 + 3.995 x 10 6 fc 4 ) _1/2 (8) 

(Ryden & Gunn 1987). The normalization constant A can be obtained 
imposing that the mass variance at 8/i~ 1 Mpc, as, is equal to unity. Filtering 
the spectrum on cluster scales, Rf = 3h~ 1 Mpc, we have obtained the rms 
torque, r(r), on a mass shell using Eq. (||) then we obtained the total 
specific angular momentum, h(r, v), acquired during expansion integrating 
the torque over time (Ryden 1988 Eq. 35): 
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Hr v) = I ( V 3 Z£^r 5/2 r {1 ~ C ° S9) m d* (9) 



where r G 5 and 5 Q are respectively the torque, the mean overdensity and 
the mean overdensity within a sphere of radius r at the current epoch to. 
The functions /bO?) are given by Ryden (1988 - Eq. 31): 

f x {Q) = 16 - 16 cos 9 + sin 2 9 -99 sin 9 (10) 

/ 2 (0) = 12- 12cos# + 3sin 2 fl-9flsin# (11) 

where 9 is a a parameter connected to the time, i, through the following 
equation: 

t = lt r o 3/ \9- sin 9) (12) 



The mean overdensity within a sphere of radius r , S(r), is given by: 



S( r , v) = — jT dxx 2 5(x) (13) 

In fig. 1 we show the variation of h(r, v) with the distance r for three values 
of the peak height v. The rms specific angular momentum, h(r, u), increases 
with distance r while peaks of greater v acquire less angular momentum 
via tidal torques. This is the angular momentum-density anticorrelation 




Figure 1: The specific angular momentum, in units of M , Mpc and the 
Hubble time, t D , for three values of the parameter v {y — 2 dotted line, 
v — 3 solid line, v — 4 dashed line) and for Rf = Sh^Mpc. 

showed by Hoffman (1986). This effect arises because the angular momen- 
tum is proportional to the gain at turn around time, t m , which in turn is 

proportional to 8(r,v)~* oc z/" 3//2 . 

3 Non-radial motions and the density criti- 
cal threshold. 

One of the consequences of the angular momentum acquisition by a mass 
shell of a proto-cluster is the delay of the collapse of the proto-structure. 



As shown by Barrow & Silk (1981) and Szalay & Silk (1983) the gravita- 
tional interaction of the irregular mass distribution of proto-cluster with 
the neighbouring proto-structures gives rise to non-radial motions, within 
the protocluster, which are expected to slow the rate of growth of the 
density contrast and to delay or suppress collapse. According to Davis & 
Peebles (1977) the kinetic energy of the resulting non-radial motions at the 
epoch of maximum expansion increases so much to oppose the recollapse of 
the proto-structure. Numerical N-body simulations by Villumsen & Davis 
(1986) showed a tendency to reproduce this so called previrialization ef- 
fect. In a more recent paper by Peebles (1990) the slowing of the growth 
of density fluctuations and the collapse suppression after the epoch of the 
maximum expansion were re-obtained using a numerical action method. In 
the central regions of a density peak (r < 0.5i?/) the velocity dispersion 
attain nearly the same value (Antonuccio & Colafrancesco 1997) while at 
larger radii (r > Rf) the radial component is lower than the tangential 
component. This means that motions in the outer regions are predom- 
inantly non-radial and in these regions the fate of the infalling material 
could be influenced by the amount of tangential velocity relative to the ra- 
dial one. This can be shown writing the equation of motion of a spherically 
symmetric mass distribution with density n(r): 

^ n{Vr) + |. n (^) + (2(vl) - (vl)) - r + n(r)^{v r ) = (14) 

where (v r ) and (v are, respectively, the mean radial and tangential stream- 
ing velocity. Eq. ( |T4"D shows that high tangential velocity dispersion 

(i v l) > 2( v r)) ma y alter the infall pattern. The expected delay in the 
collapse of a perturbation may be calculated using a model due to Peebles 
(Peebles 1993). 

We consider an ensemble of gravitationally growing mass concentrations, 
we suppose that the material in each system collects within the same po- 
tential well with inward pointing acceleration given by g{r,t). We indicate 
with dP = f(L, rv r , t)dLdv r dr the probability that a particle can be found 
in the proper radius range r, r + dr, in the radial velocity range v r = r, 
v r + dv r and with angular momentum L = tvq in the range dL. The radial 
acceleration of the particle is: 

dv r L 2 (r,u) 

where g(r) is the acceleration. Eq. ( ]T5|) can be derived from a potential 
and then from Liouville's theorem follows that the distribution function, 
/, satisfies the collisionless Boltzmann equation: 



dt r dr dv r 
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(16) 



Using Gunn & Gott's (1972) notation we write the proper radius of a shell 
in terms of the expansion parameter, a(rj,t), where is the initial radius: 



(17) 



and remembering that M = ^p(ri,t)a 3 (ri,t)rf, that |^ = p«, where p C i 
and ifj are respectively the critical mass density and the Hubble constant 
at the time ti, and assuming that no shell crossing occurs so that the total 

mass inside each shell remains constant, (p(ri,t) = ^3^!'^ ) Eq. flU 
be written as: 

d 2 a H?(l + S) AG 2 L 2 
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pi ~ pc% or integrating the equation once more: 
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where C is the binding energy of the shell. Integrating once more we have: 

iax da 
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Using Eqs (^) and ( PU| ) it is possible to find the linear over-density at 
the turn-around epoch, t ta . In fact solving Eq. (J2H) , for some epoch of 
interest, we may obtain the expansion parameter of the turn- around epoch. 
This is related to the binding energy of the shell containing mass M by 
Eq. ( |19D with ^ = 0. In turn the binding energy of a growing mode 
solution is uniquely given by the linear overdensity, 5{, at time ti. From 
this overdensity, using linear theory, we may obtain that of the turn-around 
epoch. We find the binding energy of the shell, C, using the relation 
between v and Si for the growing mode (Peebles 1980) in Eq. (|T9|) and 
finally the linear overdensity at the time of collapse: 

SG 2 L 2 ■ da] 

1 + 



S„ 



v) 



i°s(i + s) 21 1 (21) 

where S co = 1.68 is the critical threshold for a spherical model, while H 
and Q are respectively the Hubble constant and the density parameter at 
the current epoch i - We may find the angular momentum, L, needed to 
calculate S c using Eq.@. 

The mass dependence of the threshold parameter, S c (u), can be found 
as follows: we calculate the binding radius, r&, of the shell using Hoffmann 
& Shaham's criterion (1985): 

T c (r, v) < t (22) 

where T c (r, v) is the calculated time of collapse of a shell and t D is the 
Hubble time. We find a relation between v and M through the equation 
M = 47r/Vk/3. We so obtain S c (v(M)). In Fig. 2 we show the vari- 
ation of the threshold parameter, S C (M), with the mass M. Non-radial 
motions influence the value of S c increasing its value for peaks of low mass 
while leaving its value unchanged for high mass peaks. As a consequence, 
the structure formation by low mass peaks is inhibited. In other words, 
in agreement with the cooperative galaxy formation theory (Bower et al. 
1993), structures form more easily in over-populated regions. 
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Figure 2: The threshold 5 C in function of the mass M, for a CDM spectrum 
(Qo = 1, h = 1/2) with Rf = Sh^Mpc, taking account of non-radial 
motions. 

4 Tidal field and the selection function 

According to biased galaxy formation theory the sites of formation of struc- 
tures of mass ~ M must be identified with the maxima of the density peak 
smoothed over a scale R f (M oc Rj). A necessary condition for a pertur- 
bation to form a structure is that it goes nonlinear and that the linearly 
extrapolated density contrast reaches the value S(r) > 5 C = 1.68 or equiva- 
lently that the threshold criterion v t > 8 c /a (Rf) is satisfied, being a (Rf) 
the variance of the density field smoothed on scale Rf. When these con- 
dition are satisfied the matter in a shell around a peak falls in toward the 
cluster center and virializes. In this scenario only rare high v peaks form 
bright objects while low v peaks [y 1) form under-luminous objects. 
The kind of objects that form from nonlinear structures depends on the 
details of the collapse. Moreover if structures form only at peaks in the 
mass distribution they will be more strongly clustered than the mass. Sev- 
eral feedback mechanisms has been proposed to explain this segregation 
effect (Rees 1985, Dekel & Rees 1987). Even if these feedback mechanisms 
work one cannot expect they have effect instantaneously, so the threshold 
for structure formation cannot be sharp (BBKS). To take into account this 



effect BBKS introduced a threshold or selection function, t{y/v t ). The se- 
lection function, t(i//i/ t ), gives the probability that a density peak forms 
an object, while the threshold level, u t , is defined so that the probability 
that a peak form an object is 1/2 when v = v t . The selection function 
introduced by BBKS (Eq. 4.13), is an empirical one and depends on two 
parameters: the threshold v t and the shape parameter q: 

1 + {vjvtf 

If q — > oo this selection function is a Heaviside function d(v — u t ) so that 
peaks with v > v t have a probability equal to 100% to form objects while 
peaks with v < v t do not form objects. If q has a finite value sub-z/ 4 peaks 
are selected with non-zero probability. Using the given selection function 
the cumulative number density of peaks higher than v is given, according 
to BBKS, by: 



n pk = / t(is/is t )N pk (v)dv (24) 

where N p k{v) is the comoving peak density (see BBKS Eq. 4.3). A form of 
the selection function, physically motivated, can be obtained following the 

argument given in CAD. In this last paper the selection function is defined 
as: 



t{v) = I p 5,{5)(r M t,v),Vs( r Mt,v) d5 (25) 
where the function 
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gives the probability that the peak overdensity is different from the average, 
in a Gaussian density field. The selection function depends on v through 
the dependence of 5(r) from v. As displayed the integrand is evaluated at 
a radius rut which is the typical radius of the object that we are selecting. 
Moreover the selection function tiv) depends on the critical overdensity 
threshold for the collapse, 5 C , which is not constant as in a spherical model 
(due to the presence, in our analysis, of non-radial motions that delay the 
collapse of the proto-cluster) but it depends on v. Known 5 c {v) and chosen 
a spectrum, the selection function is immediately obtainable through Eq. 
(|25| ) and Eq. (p6|). The result of the calculation, plotted in Fig. 3, for 
two values of the filtering radius, (Rf = 2, 3 h~ l Mpc), shows that the 
selection function, as expected, differs from an Heaviside function (sharp 
threshold). The value of v at which the selection function tiv) reaches the 
value 1 {t(v) ~ 1) increases for growing values of the filtering radius, Rf. 
This is due to the smoothing effect of the filtering process. The effect of 
non-radial motions is, firstly, that of shifting t(z/) towards higher values 
of u, and, secondly, that of making it steeper. The selection function is 
also different from that used by BBKS (Tab. 3a). Finally it is interesting 
to note that the selection function defined by Eq. (|25|) and Eq. fl26|) is 
totally general, it does not depend on the presence or absence of non- 
radial motions. The latter influence the selection function form through 
the changement of 5 C induced by non-radial motions itself. 




Figure 3: The selection function, t(u), for Rf = 3h~ 1 Mpc (S c = 1.68, solid 
line; S c function of u, dotted line) and for 4h~ 1 Mpc (5 C = 1.68, short dashed 
line; 5 C function of v, long dashed line). 

5 The bias coefficient 

A model of the Universe in which light traces the mass distribution accu- 
rately (unbiased model) is subject to several problems. As pointed out by 
Davis et al. (1985) an unbiased CDM produces a galaxy correlation func- 
tion which is steeper than observed and a pairwise velocity dispersion larger 
than that deduced from redshift surveys. A remedy to this problem can 
be found if the assumption that light trace mass is relaxed introducing the 
biasing concept, i.e. that galaxies are more clustered than the distribution 
of matter in agreement to the concept of biasing inspired by Kaiser's (1984) 
discussion of the observation that clusters of galaxies cluster more strongly 
than do galaxies, in the sense that the dimensionless two-point correlation 
function, £cc(r), is much larger than the galaxy two-point function, £ gg (r). 
The galaxy two-point correlation function £ gg (r) is a power-law : 

Ur) = (— ) 7 (27) 

r 0,g 

with a correlation length r 0)9 ~ 5/i~ x Mpc and a slope 7 ~ 1.8 for r < 10/z 1 
Mpc (Davis & Peebles 1983; Davis et al. 1985; Shanks et al. 1989), (some 



authors disagree with this values; for example Strauss et al. 1992 and 
Fisher et al. 1993 find ro )S — 3.79/iT 1 Mpc and 7 ~ 1.57). As regards the 
clusters of galaxies the form of the two-point correlation function, £ cc (r), is 
equal to that given by Eq. (p7|). Only the correlation length is different. In 
the case of clusters of galaxies the value of ro, c is uncertain (see Bahcall & 
Soneira 1983; Postman et al. 1986; Sutherland 1988; Bahcall 1988; Dekel 
et al. 1989; Olivier et al. 1990 and Sutherland & Efstathiou 1991) however 
it lays in the range ro, c — 12-^25/i _1 Mpc in any case larger than ro i9 . One 
way of defining the bias coefficient of a class of objects is that given by 
(BBKS): 

b(R f ) = M + 1 (28) 



where (z/) is: 
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t(-)N pk (u)du (29) 



from Eq. ( p9"D it is clear that the bias parameter can be calculated once 
a spectrum, P(k), is fixed. The bias parameter depends on the shape 
and normalization of the power spectrum. A larger value is obtained for 
spectra with more power on large scale (Kauffmann et al. 1996). In this 
calculation we continue to use the standard CDM spectrum (fl Q = 1, h = 
0.5) normalized imposing that the rms density fluctuations in a sphere 
of radius 8h~ 1 Mpc is the same as that observed in galaxy counts, i.e. 
erg = a(8h~ 1 Mpc) = 1. The calculations have been performed for three 
different values of the filtering radius (Rf = 2, 3, 4 h~ l Mpc). The values 
of b, that we have obteined, are respectively, in order growing of Rf, 1.6, 
1.93 and 2.25. 

As shown, the value of the bias parameter tends to increase with Rf 
due the filter effect of tiy). As shown t(u) acts as a filter, increasing the 
filtering radius, Rf, the value of v at which t(v) ~ 1 increases . In other 
words when Rf increases t(v) selects density peaks of larger height. The 
reason of this behavior must be searched in the smoothing effect that the 
increasing of the filtering radius produces on density peaks. When Rf is 
increased the density field is smoothed and t(z/) has to shift towards higher 
value of v in order to select a class of object of fixed mass, M. 



6 The X-ray temperature function 

The PS theory provides an analytical description of the evolution of struc- 
tures in a hierarchical Universe. In this model the linear density field, 
p(x, t), is an isotropic random Gaussian field, the non-linear clumps are 
identified as over-densities (having a density contrast 5 C ~ 1.68 - Gunn 
& Gott 1972) in the linear density field, while a mass element is incorpo- 
rated into a non-linear object of mass M when the density field smoothed 
with a top-hat filter of radius Rf, exceeds a threshold 5 C (M oc R 3 f). The 



probability distribution for fluctuations is given by: 



p[6(M)} = exp[-5(M) 2 /2a(M) 2 ] (30) 

The probability that an object of mass M has formed is obtained inte- 
grating Eq. (|30D from the threshold value 5 C to infinity and the comoving 
number density of non-linear objects of mass M to M + dM is given simply 
by differentiating the integral with respect to mass and is given by: 

N(M, tyi M = -„^exp (-.72) \ (£) f (31) 

where pb is the mean mass density, a(M) is the rms linear mass overdensity 
evaluated at the epoch when the mass function is desidered and v = -^f^- 
The redshift dependence of Eq. ( |3~T1) can be obtained remembering that 



_ 5 c (z)D(0) 

~ a (M)D(z) {62} 

being D(z) the growth factor of the density perturbation and <J (M) the 
current value of a(M). In Eq. ([31]) PS introduced arbitrarily a factor of 
two because J °° dF(M) = 1/2, so that only half of the mass in the Universe 
is accounted for. Bond et al. (1991) showed that the "fudge factor" 2 is 
naturally obtained using the excursion set formalism in the sharp fc-space 
while for general filters (e.g., Gaussian or "top hat") it is not possible to 
obtain an analogous analytical result. As stressed by Yano et al. 1996, the 
factor of 2 obtained in the sharp fc-space is correct only if the spatial corre- 
lation of the density fluctuations is neglected. In spite of the quoted prob- 
lem, several authors (Efstathiou et al. 1988; Brainerd & Villumsen 1992; 
Lacey & Cole 1994) showed that PS analytic theory correctly agrees with 
N-body simulations. In particular Efstathiou et al. (1988), showed that 
PS theory correctly agrees with the evolution of the distribution of mass 
amongst groups and clusters of galaxies (multiplicity function). Brainerd 
& Villumsen (1992) studied the CDM halo mass function using a hierar- 
chical particle mesh code. From this last work results that PS formula fits 
the results of the simulation up to a mass of 10 times the characteristic 
la fluctuation mass, M*, being M* ~ lO 15 &~ 6 /( ni+3 )M , where b is the bias 
parameter and ni is the local slope of the power spectrum. PS theory has 
proven particularly useful in analyzing the number counts and redshift dis- 
tributions for QSOs (Efstathiou & Rees 1988), Lyman a clouds (Bond et 
al. 1988) and X-ray clusters (Cavaliere & Colafrancesco 1988). 
Some difficulties arise when PS theory is compared with observed distri- 
butions. To estimate the multiplicity function of real systems it is in fact 
required a knowledge of the temperature-mass (T-M) relation in order to 
trasform the mass distribution into the temperature distribution. Theo- 
retical uncertainty arises in this transformation because the exact relation 
between the mass appearing in the PS expression and the temperature of 
the intracluster gas is unknown. Under the standard assumption of the 
Intra-Cluster (IC) gas in hydrostatic equilibrium with the potential well 



of a spherically simmetric, virialized cluster, the IC gas temperature-mass 
relation is easily obtained by applying the virial theorem and for a flat 
matter-dominated Universe we have that (Kaiser 1986, Evrard 1990): 



T = (B.4h><*keV) (j^- J (1 + ,) (33) 

The assumptions of perfect hydrostatic equilibrium and virialization are 
in reality non completely satisfied by clusters. Clusters profile may have 
departure from isothermality, with sligth temperature gradients throughout 
the cluster. The X-ray weighted temperature can be slightly different from 
the mean mass weighted virial temperature. In any case the scatter in the 
T-M relation given by Eq. (^) is of the order of ~ 10% (Evrard 1991). As 
shown by Bartlett & Silk (1993) the X-ray distribution function obtained 
using a standard CDM spectrum over-produces the clusters abundances 
data obtained from Henry & Arnaud (1991) and Edge et al. (1990). The 
discrepancy can be reduced taking into account the non-radial motions 
that originate when a cluster reaches the non-linear regime. In fact, the PS 
temperature distribution requires specification of 5 C and the temperature- 
mass relation T-M. The presence of non-radial motions changes both 5 C 
and the T-M relation. To get the temperature distribution it is necessary 
to know the temperature-mass relation. This can be obtained using the 
virial theorem, energy conservation and using Eq. ( |i~9"|) (Bartlett & Silk 
1993). From the virial theorem we may write: 

, v GM reff L 2 

(K) = + / — dr 34 

V 1 2r eff Jo 2M 2 r 3 v ; 

while from the energy conservation: 

tT ^ GM reff L 2 , GM rn« L 2 , 
-(K) + + / — - -dr = + / -—dr 35 

r e ff Jo M z r 6 r ta Jo M z r- i 

Eq. (|3~4D and Eq.([35|) can be solved for r e ff and < K >. We finally have 
that: 



T = (6AkeV) 



M-h \ 2/3 



1O 15 M 



I i Jo M2% 

( G 2^£ M 2)l/3 



(36) 



where rj is a parameter given by rj = r ta /x\, being r ta the radius of the 
turn-around epoch, while x± is defined by the relation M = Airpi ) x\/?> and 
if; = r e ff/r ta where r e ff is the time-averaged radius of a mass shell. Eq. 
(|36| ) was normalised to agree with Evrard's (1990) simulations for L = 0. 
The new T-M relation, Eq. (|36|), differs from Eq. (|33| ) for the presence of 



the term: 



Jo M^r 3 



(37) 



This last term changes the dependence of the temperature from the mass, 
M, in the T-M relation. Moreover the new T-M relation depends on the 



angular momentum, L, originating from the gravitational interaction of the 
quadrupole moment of the protocluster with the tidal field of the matter 
of the neighboring protostructures. In Fig. 4 the X-ray temperature distri- 
bution, derived using a CDM model with Q Q = 1, h = 1/2 and taking into 
account of non-radial motions, is compared with Henry & Arnaud (1991) 
and Edge et al. (1990) data and with a pure CDM model with Q = 1, 
h — 1/2. As shown the CDM model that does not take account of the non- 




Figure 4: X-ray temperature distribution function. The solid line gives 
the temperature function for a pure CDM model (Qq = 1, /i = 1/2), with 
Rf = ?>hr x M'pc. The dotted line is the same distribution but now taking 
account of non-radial motions. The data are obtained by Edge et al. 1990, 
and Henry & Arnaud 1991 

radial motions over-produces the clusters abundance. The introduction of 
non-radial motions gives a more careful description of the experimental 
data. As we have seen the X-ray temperature distribution function ob- 
tained taking account of non-radial motions is different from that of a pure 
CDM model for two reasons: 

1) the variation of the threshold, <5 C , with mass, M. This is due to the 
changement of the energetics of the collapse model produced by the intro- 
duction of another potential energy term { ^fitI ) in -Eg. ([15]); 

2) the modification of the T-M relation produced by the alteration of the 
partition of energy in virial equilibrium. 



For values of mass M = O.5M the difference between the two theoretical 
lines in Fig. 4 is due to the first factor for a ~ 59% and this value increases 
with increasing mass. The uncertainty in our model fundamentally comes 
from the uncertainty of the T-M relation whose value has been previously 
quoted. 

One of the objection to the result of this paper may that the effect de- 
scribed has not been seen in some hydrodynamic simulations (Evrard & 
Crone 1992). The answer to this objection is that our model is funda- 
mentally based on the previrialization conjecture (Davis & Peebles 1977; 
Peebles 1990), (supposing that initial asphericities and tidal interactions 
between neighboring density fluctuations induce significant non-radial mo- 
tions, which oppose the collapse) and while some N-body simulations (Vil- 
lumsen & Davis 1986; Peebles 1990) appear to reproduce this effect, other 
simulations (for example Evrard & Crone 1992) do not. An answer to this 
controversy was given by (Lokas et al. 1996). The problem is connected to 
spectral index n used in the simulations. The "previrialization" is seen only 
for n > — 1. While Peebles (1990) used simulations with n = 0, Evrard & 
Crone (1992) assumed n = — 1. Excluding this particular case generally 
the ensemble properties of clusters like their optical and X-ray luminosity 
functions, or their velocity and temperature distribution functions, are dif- 
ficult to address directly in numerical simulations because the size of the 
box must be very large in order to contain a sufficient number of clusters; 
an then analytical approach remains an effective alternative. 



7 Non-radial motions and the clusters cor- 
relation function 

Supposing that the Universe at some early epoch can be described using 
the same assumptions of the previous section the probability p(M\, M 2 ,r) 
per unit volume per unit masses of finding two collapsed objects of mass M\ 
and M 2 separated by a distance r is obtained by integrating both variables 
of the bivariant Gaussian distribution in 8(Mi) and 8{M2), with correlation 
£ p (r) from 5 C to infinity and then taking the partial derivatives with respect 
to both masses. The correlation function for collapsed objects is simply 
obtained from: 

Cmm = p{Mt, M 2 , r)/p(M 1 )p(M 2 ) - 1 (38) 
which for equal masses and weak correlations £ p <C 1 is (Kashlinsky 1987): 

tMM = [5 2 Ja p (M) 2 ]Ur) (39) 

where a p (M) is the variance of the mass fluctuation and £ p is the correlation 
function of the matter density distribution when the density fluctuations 
had small amplitude. Eq. (^) shows that the correlation of collapsed ob- 
jects may be enhanced relative to that of the underlying mass fluctuations. 
This condition is usually described by the bias parameter b which is some- 
times defined as [CMM(r)/^(r)] 1/2 = 5 c /a p (M). 

Studies of clustering on scales > 10/i -1 Mpc have shown that the correlation 



function given in Eq. ( p9|) is different from that obtained from observa- 
tions. The most compelling data are angular correlation functions for the 
APM survey. These decline much less rapidly on large scales than the CDM 
prediction (Maddox et al. 1990). As discussed in the introduction there 
are two ways to reduce the quoted discrepance: either with a modification 
of CDM theory involving the physics of early Universe, or with a modi- 
fication of CDM theory involving the physics of galaxy formation. This 
discrepancy can be reduced, similarly to the problem of over- abundance of 
clusters, taking into account the non-radial motions that originate when a 
cluster reaches the non-linear regime. In fact, the calculation of the correla- 
tion function requires the specification of S c which is changed by non-radial 
motions and whose new value is given by Eq. (pip. In Fig. 5 the correla- 
tion function derived taking into account non-radial motions is compared 
both with the two-point correlation function obtained by Sutherland & Ef- 
stathiou (1991) from the analysis of Huchra's et al. (1990) deep redshift 
survey as discussed in Geller & Huchra (1988) and with the data points 
for the APM clusters computed by Efstathiou et al. (1992). As shown 
the discrepancy between pure CDM previsions and experimental data is 
remarkable. The CDM model seems to have trouble in re-producing the 
behaviour of the data. In fact, the predicted two-point cluster function 
is too steep and rapidly goes nearly to zero for r ~ 30h~ 1 Mpc, while the 
data show no significant anticorrelation up to r ~ 60/i -1 Mpc (see Borgani 
1990). The introduction of non-radial motions gives a more accurate de- 
scription of the experimental data. The result obtained is in agreement 
with that of Borgani (1990) who studied the effect of particular thresh- 
olds (erfc-threshold and Gaussian-threshold) on the correlation properties 
of clusters of galaxies. The fundamental difference between our and Bor- 
gani's approach is that our threshold function is physically motivated: it is 
simply obtained from the assumptions of a Gaussian density field and tak- 
ing account of non-radial motions. Borgani's threshold functions (erfc and 
Gaussian threshold) are ad-hoc introduced in order to reduce the discrep- 
ancy between the observed and the CDM predicted two points correlation 
functions of clusters of galaxies. The connection with the quoted non- 
sphericity effects, even if logical and in agreement with our results, is only 
a posteriori tentative to justify the choice made. 



8 Conclusions 

In these last years many authors have shown the existence of a strong dis- 
crepancy between the observed properties of clusters of galaxies and that 
predicted by the CDM model. To reduce this discrepancy several alterna- 
tive models have been introduced but no model has considered the role of 
the non-radial motions. Here we have shown how non-radial motions may 
reduce some of these discrepancies: namely 

a) the discrepancy between the observed £(r) of cluster of galaxies and that 
predicted by CDM model; 

b) the discrepancy between the measured X-ray temperature distribution 
function and that predicted by the CDM model. 

To this aim we calculated the variation of the treshold parameter, 5 C , as 
a function of the mass M, produced by the presence of non-radial mo- 
tions in the outskirts of clusters of galaxies. We used S C (M) to calculate 




r (Mpc) 



Figure 5: Clusters of galaxies correlation function. The solid line gives 
the correlation function for a pure CDM model, with R f = ZhT^Mpc. 
The dashed line is the same distribution but now taking account of non- 
radial motions. The observational data refer to the two-point correlation 
function obtained by Sutherland & Efstathiou (1991) {filled exagons) from 
the analysis of Huchra's et al. (1990) deep redshift survey and with the 
data points for the APM clusters computed by Efstathiou et al. (1992) 
{dashed errorbars). 



the two-point correlation of clusters of galaxies and the X-ray temperature 
distribution function. We compared the prediction for the two point corre- 
lation function with that obtained by Sutherland & Efstathiou (1991) from 
the analysis of Huchra's et al. (1990) deep redshift survey as discussed in 
Geller & Huchra (1988) and with the data points for the APM clusters 
computed by Efstathiou et al. (1992). The prediction for the X-ray tem- 
perature function was compared with Henry & Arnaud (1991) and Edge et 
al. (1990) X-ray temperature distributions for local clusters. Our results 
(see Fig. 4 and Fig. 5) show how the non-radial motions change both the 
correlation length of the correlation function, making it less steep than that 
obtained from a pure CDM model where the non-radial motions are not 
considered, and the X-ray temperature function. In both cases our model 
gives a good agreement with the data. 

Finally we calculated the bias coefficient using a selection function that 
takes into account the effects of non-radial motions, and we show that the 
bias so obtained can account for a substantial part of the total bias required 
by observations on cluster scales. 
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